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ABSTRACT 

Complementary views of galaxy clusters in the radio synchrotron, hard X-ray inverse Comp- 
ton, and high-energy y-ray regimes are critical in calibrating them as high-precision cosmo- 
logical probes. We present predictions for scaling relations between cluster mass and these 
non-thermal observables. To this end, we use high-resolution simulations of a sample of 
galaxy clusters spanning a mass range of almost two orders of magnitudes, and follow self- 
consistent cosmic ray physics on top of the radiative hydrodynamics. We model relativistic 
electrons that are accelerated at cosmological structure formation shocks and those that are 
produced in hadronic interactions of cosmic rays with ambient gas protons. Calibrating the 
magnetic fields of our model with Faraday rotation measurements, the synchrotron emission 
of our relativistic electron populations matches the radio synchrotron luminosities and mor- 
phologies of observed giant radio halos and mini-halos surprisingly well. Using the complete 
sample of the brightest X-ray clusters observed by ROSAT in combination with our y-ray 
scaling relation, we predict GLAST that will detect about ten clusters allowing for Eddington 
bias due to the scatter in the scaling relation. The expected brightest y-ray clusters are Ophi- 
uchus, Fornax, Coma, A3627, Perseus, and Centaurus. The high-energy y-ray emission above 
100 MeV is dominated by pion decays resulting from hadronic cosmic ray interactions. We 
provide an absolute lower flux limit for the y-ray emission of Coma in the hadronic model 
which can be made tighter for magnetic field values derived from rotation measurements to 
match the GLAST sensitivity, providing thus a unique test for the possible hadronic origin of 
radio halos. Our predicted hard X-ray emission, due to inverse Compton emission of shock 
accelerated and hadronically produced relativistic electrons, falls short of the detections in 
Coma and Perseus by a factor of 50. This casts doubts on inverse Compton interpretation and 
reinforces the known discrepancy of magnetic field estimates from Faraday rotation measure- 
ments and those obtained by combining synchrotron and inverse Compton emission. 

Key words: galaxies: cluster: general, cosmic rays, magnetic fields, radiation mechanisms: 
non-thermal 



1 INTRODUCTION 

Previously, it has been assumed that galaxy clusters are sufficiently 
well described by their mass which was thought to be largely inde- 
pendent of the complex astrophysical processes taking place in the 
intra-cluster medium (ICM) such as star formation and different 
kinds of feedback processes. High-resolution XMM-Newton and 
Chandra X-ray observations taught us in the last years that this 
over-simplified paradigm needs to be modified. Even ostensibly 
'relaxed' clusters reveal a richness of substructure with substan- 
tial small-scale variation in temperature, metallicity, and surface 
brightness. 

This raises the question if high-precision cosmology will prin- 
cipally be possible using clusters. Clearly, we need to understand 
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how non-equilibrium processes that lead to cosmic ray populations 
and turbulence impact on the thermal X-ray emission and Sunyaev- 
Zel'dovich effect. This forces us to explore complementary obser- 
vational windows to clusters such as non-thermal emission that can 
potentially elucidate the otherwise invisible non-equilibrium pro- 
cesses. The upcoming generation of low-frequency radio, hard X- 
ray, and y-ray instruments open up the extragalactic sky in unex- 
plored wavelength ranges (cf. Pfrommer et al. 2007b, for a compi- 
lation of these experiments). Suitably combining radio synchrotron 
radiation, inverse Compton emission in the hard X-ray regime, and 
high-energy y-ray emission will enable us to estimate the cosmic 
ray pressure contribution and provide us with clues to the dynami- 
cal state of a cluster. This will allow us to construct a 'gold sample' 
for cosmology using information on the dynamical cluster activ- 
ity that is orthogonal to the thermal cluster observables. Addition- 
ally, these non-thermal observations have the potential to improve 
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our knowledge of diffusive shock acceleration, large scale magnetic 
fields, and turbulence. 

Of the possible non-thermal emission bands, only the dif- 
fuse large-scale radio synchrotron emission of clusters has been 
unambiguously detected so far. Generally these radio phenomena 
can be divided into two categories that differ morphologically, in 
their degree of polarisation, as well as in their characteristic emis- 
sion regions with respect to the cluster halo. The large-scale 'ra- 
dio relic' or 'radio gischt' emission (Kempner et al. 2004), that 
has a high degree of polarisation, is irregularly shaped and occurs 
at peripheral cluster regions, can be attributed to merging or ac- 
cretion shock waves as proposed by Ensslin et al. (1998). Promi- 
nent examples for large scale 'radio relic' emission have been ob- 
served in Abell 3667 (Rottgering et al. 1997), Abell 3376 (Bagchi 
et al. 2006), and Abell 2256 (Clarke & EnlSlin 2006). In contrast, 
'cluster radio halos' show a coherently large diffuse radio emis- 
sion that is centred on the cluster, resemble the underlying ther- 
mal bremsstrahlung emission in X-rays, are unpolarised, and show 
spectral index variations that are amplified in the peripheral regions 
of the extended radio emitting regions. These radio halo phenom- 
ena can be furthermore subdivided into Mpc-sized 'giant radio ba- 
les' that are associated with merging clusters and 'radio mini-halos' 
that are observed in a few cool core clusters and have a smaller ex- 
tent of a few hundreds of kpc. Prominent examples for 'giant radio 
halos' can be obtained from Giovannini et al. (1999) and include 
the Coma cluster (Kim et al. 1989; Deiss et al. 1997) and the galaxy 
cluster IE 0657-56 (Liang et al. 2000). Prominent 'radio mini ha- 
los' are observed in the Perseus cluster (Pedlar et al. 1990) or in 
RX J1347.5-1 145 (Gitti et al. 2007). 

Previously, there have been two models suggested that are able 
to explain 'cluster radio halos'. (1) Re -acceleration processes of 
'mildly' relativistic electrons (y - 100 - 300) that are being in- 
jected over cosmological timescales into the ICM by sources like 
radio galaxies, merger shocks, or galactic winds can provide an ef- 
ficient supply of highly-energetic relativistic electrons. Owing to 
their long lifetimes of a few times 10' years these 'mildly' rela- 
tivistic electrons can accumulate within the ICM (Sarazin 2002), 
until they experience continuous in- situ acceleration either via in- 
teractions with magneto-hydrodynamic waves, or through turbulent 
spectra (Jaffe 1977; Schlickeiser et al. 1987; Brunetti et al. 2001; 
Ohno et al. 2002; Brunetti et al. 2004; Gitti et al. 2004; Brunetti 
& Lazarian 2007). (2) Hadronic interactions of relativistic protons 
with ambient gas protons produce pions which decay successively 
into secondary electrons, neutrinos and y-rays. These secondary 
relativistic electrons and positrons can emit a halo of radio syn- 
chrotron emission in the presence of ubiquitous intra-cluster mag- 
netic fields (Dennison 1980; Vesttand 1982; Blasi & Colafrancesco 
1999; Dolag & EnBlin 2000; Miniati et al. 2001b; Pfrommer & 
EnBlin 2003, 2004a,b) as well as inverse Compton emission by 
scattering photons from the cosmic microwave background into the 
hard X-ray and y-regime. In our companion paper (Pfrommer et al. 
2007b, hereafter Paper II), we suggest a modification of the lat- 
ter model that is motivated by our high-resolution cluster simu- 
lations and cures the weaknesses of the original model. We find, 
that our simulated giant radio halos are dominated in the centre 
by secondary synchrotron emission with a ttansition to the radio 
synchrotton radiation emitted from shock-accelerated electrons in 
the cluster periphery. This explains the extended radio emission 
found in merging clusters, while it is more centrally concentrated 
in relaxed cool core clusters. Varying spectral index distributions 
preferably in the cluster periphery (Feretti et al. 2004) support this 
picture. The characterisation of quantities related to cosmic rays in 



clusters can be found in our first companion paper that studies the 
interplay of thermal gas and cosmic rays and their effect on ther- 
mal cluster observables such as X-ray emission and the Sunyaev- 
Zel'dovich effect (Pfrommer et al. 2007a, hereafter Paper I). 

The outline of the paper is as follows. Section 2 describes our 
general methodology, presents our cluster sample, and the different 
simulated physical processes. In Sect. 3, we present the results on 
the cluster scaling relations for non-thermal observables as well as 
y-ray flux and luminosity functions. These are compared to obser- 
vations and finally critically discussed in Sect. 4. 



2 METHODOLOGY 
2.1 General procedure 

We have performed high-resolution hydrodynamic simulations of 
the formation of 14 galaxy clusters. The clusters span a mass range 
from 5 X 10'^ ft ' Mq to 2 x 10'^ ft"' Mq and show a variety of dy- 
namical states ranging from relaxed cool core clusters to violent 
merging clusters (cf. Table 1). Our simulations dynamically evolve 
dissipative gas physics including radiative cooling, star formation, 
and supernova feedback. We identify the strength of structure for- 
mation shock waves on-the-fly in our simulations and measure the 
shock Mach number that is defined by the ratio of shock velocity 
to pre-shock sound velocity, M. = f'shock/Csound (Pfrommer et al. 
2006). On top of this, we self-consistently follow cosmic ray (CR) 
physics including adiabatic CR transport processes, injection by su- 
pemovae and cosmological structure formation shocks, as well as 
CR thermalization by Coulomb interaction and catastrophic losses 
by hadronic interactions (EnlSlin et al. 2007; Jubelgas et al. 2007). 
In our post-processing, we model relativistic electrons that are ac- 
celerated at cosmological structure formation shocks and those that 
are produced in hadronic interactions of cosmic rays with ambi- 
ent gas protons. This approach is justified since these elections do 
not modify the hydrodynamics of the gas owing to their negligible 
pressure contribution. We compute the stationary relativistic elec- 
tron spectrum that is obtained by balancing the mentioned injection 
mechanisms with the synchrotton and inverse Compton cooling 
processes. Details of our modelling can be fovmd in Paper n. Both 
populations of relativistic electrons emit a morphologically distin- 
guishable radio synchrotton radiation as well as inverse Compton 
emission due to up-scattering of photons of the cosmic microwave 
background (CMB) into the hard X-ray and y-ray regime. At en- 
ergies larger than 100 MeV, we expect additionally y-ray emission 
from decaying pions that are produced in hadronic CR interactions. 
While the emission of the shock accelerated primary electrons is 
amorphous and peripheral as observed in radio relics, the hadron- 
ically produced secondary electrons show a centrally concentrated 
emission characteristic that resembles that of the central parts of 
observed radio halos. 

In this paper, we concenttate on three observationally moti- 
vated wave-bands. (1) Radio ssmchrotton emission at 1.4 GHz, (2) 
non-thermal hard X-ray emission at energies Ey > \Q keV, and 
(3) y-ray emission at energies Ey > 100 MeV. We study the con- 
tribution of the different emission components to the total cluster 
luminosity in each of these bands, derive cluster scaling relations, 
and study their dependence on the simulated physics and adopted 
parametrisation of the magnetic field. The radio synchrotron scal- 
ing relation is then compared to the observed sample of giant radio 
halos and radio mini-halos. Using cluster masses from the complete 
sample of the X-ray brightest clusters (HIFLUGCS, tiie Highest X- 
ray FLUx Galaxy Cluster Sample, Reiprich & Bohringer (2002)), 
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Table 1. : Cluster sample 



Cluster 


sim.'s 


1 ..in 
dyn. state^ ' 


M 

[r^ Mo] 


d(2) 
^200 

[h~^ Mpc] 


'^^ 200 

[keV] 


1 


28a 


CC 


1.8 X lO'^ 


2.0 


13.1 


2 


ffl a 


cc 


1.3 X lO'^ 


1.8 


10.6 


3 


g72a 


PostM 


1.1 X lO'^'' 


1.7 


9.4 


4 


g51 


CC 


1.1 X lO'' 


1.7 


9.4 


5 


elb 


M 


3.7 X lO'** 


1.2 


4.7 


6 


g72b 


M 


1.5 X 10" 


0.87 


2.4 


7 


glc 


M 


1.4 X 10" 


0.84 


2.3 


8 


g8b 


M 


1.0 X 10" 


0.76 


1.9 


9 


gld 


M 


9.2 X 10" 


0.73 


1.7 


10 


g676 


CC 


8.8 X 10" 


0.72 


1.7 


11 


g914 


CC 


8.5 X 10" 


0.71 


1.6 


12 


gle 


M 


6.4 X 10" 


0.65 


1.3 


13 


gSc 


M 


5.9 X 10" 


0.63 


1.3 


14 


gSd 


PreM 


5.4 X 10" 


0.61 


1.2 


Notes: 



(1) The dynamical state has been classified through a combined criterion 
invoking a merger tree study and the visual inspection of the X-ray bright- 
ness maps. The labels for the clusters are M-merger, PostM-post merger 
(slightly elongated X-ray contours, weak cool core region developing), 
PreM-pre-merger (sub-cluster already within the virial radius), CC-cool 
core cluster with extended cooling region (smooth X-ray profile). 

(2) The virial mass and radius are related by Ma(z) = Ap,.,-i,(z)R^, 
where A = 200 denotes a multiple of the critical overdensity Paitd) = 
3H(zf/(SnG). 

(3) The virial temperature is defined by kT/^ = GM/^jxt7i^/(2R/^), where 
denotes the mean molecular weight. 

we construct luminosity and flux functions for the hard X-ray and 
y-ray band. This allows us to identify the brightest clusters in the 
hard X-ray and y-ray slsy and predict the cluster sample to be seen 
by GLAST. 

2.2 Adopted cosmology and cluster sample 

We provide only a short overview over the simulations and our 
cluster sample for completeness while the simulation details can be 
found in Paper II. Simulations were performed using the 'concor- 
dance' cosmological cold dark matter model with a cosmological 
constant (ACDM). The cosmological parameters of our model are: 
= ^DM + fib = 0.3, Qb = 0.039, a^ = 0.7, h = 0.7, n = 1, and 
(Tg = 0.9. Here, Q.^ denotes the total matter density in units of the 
critical density for geometrical closure today, pcru = 3Hg/(SnG). 
Ob and CIa denote the densities of baryons and the cosmological 
constant at the present day. The Hubble constant at the present day 
is parametrised as Hg = 100 h km s^'Mpc"', while n denotes the 
spectral index of the primordial power-spectrum, and iTg is the rms 
linear mass fluctuation within a sphere of radius 8 A" 'Mpc extrapo- 
lated to z = 0. 

We analysed the clusters with a halo-finder based on spherical 
overdensity followed by a merger tree analysis in order to get the 
mass accretion history of the main progenitor. We also produced 
projections of the X-ray emissivity at redshift z = in order to 
get a visual impression of the cluster morphology. The dynamical 
state of a cluster is defined by a combined criterion: (i) if the clus- 
ter did not experience a major merger with a progenitor mass ratio 
1:3 or larger after z = 0.8 (corresponding to a look-back time of 
~ 5 /i ' Gyr) and (ii) if the visual impression of the cluster's X-ray 
morphology is relaxed, it was defined to be a cool core cluster. The 



Table 2. ; Different physical processes m our SMULmoNs: 



Simulated physics*^'^ 


simulation models' 
,SI S2 


1). 

.S3 


thermal shock heating 






■/ 


radiative cooling 






•/ 


star formation 




/ 


•/ 


Coulomb CR losses 




~r 


•/ 


hadronic CR losses 


/ 


/ 




shock-CRs 


/ 






supemova-CRs 









Notes: 

(1) This table serves as an overview over our simulated models. The first 
column shows the simulated physics and the following three columns show 
our difli'erent simulation models with varying gas and cosmic ray physics. 
Model SI models the thermal gas non-radiatively and includes CR physics, 
while the models S2 and S3 use radiative gas physics with different variants 
of CR physics. 



spherical overdensity definition of the virial mass of the cluster is 
given by the material lying within a sphere centred on a local den- 
sity maximum, whose radial extend i?A is defined by the enclosed 
threshold density condition M(< i?A)/(4:T/?^/3) = pthres- We chose 
the threshold density Pihres(z) = Apcrit(z) to be a multiple A = 200 
of the critical density of the universe Pcritfe) = 3//(z)^/(8:tG). We 
assume a constant A = 200 although some treatments employ a 
time-varying A in cosmologies with \ (Eke et al. 1996). In 

the reminder of the paper, we use the terminology i?vii instead of 

^200- 



2.3 The models 

For each galaxy cluster we ran three different simulations with 
varying gas and cosmic ray physics (cf. Table 2). The first set of 
simulations used non-radiative gas physics only, i.e. the gas is trans- 
ported adiabatically unless it experiences structure formation shock 
waves that supply the gas with entropy and thermal pressure sup- 
port. Additionally we follow cosmic ray (CR) physics including 
adiabatic CR transport processes, injection by cosmological struc- 
ture formation shocks with a Mach number dependent acceleration 
scheme, as well as CR thermalization by Coulomb interaction and 
catastrophic losses by hadronic interactions (model SI). The sec- 
ond set of simulations follows the radiative cooling of the gas, star 
formation, supernova feedback, and a photo-ionising background 
(details can be found in Paper I). As before in model S 1 . wc ac- 
count for CR acceleration at structure formation shocks and allow 
for all CR loss processes (model S2). The last set of simulations ad- 
ditionally assumes that a constant fraction = «^CR,mj/£diss = 0.3 
of the kinetic energy of a supernova ends up in the CR population 
(model S3), which is motivated by TeV y-ray observations of a su- 
pernova remnant that find an energy fraction of ^sn - 0.1 - 0.3 
when extrapolating the CR distribution function (Aharonian et al. 
2006). We choose a maximum value for this supernova energy ef- 
ficiency owing to the large uncertainties and our aim to bracket the 
realistic case with the two radiative CR simulations. 

Since we have not run self-consistent magneto-hydrodynam- 
ical (MHD) simulations on top of the radiative gas and CR physics, 
we chose the following model for the magnetic energy density to 
compute the synchrotron and inverse Compton (IC) emission: 
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I, Eth.O 



Tan 



(1) 



Synchrotron emission (v — 1.4 GHz): 



where the central magnetic energy density sb,o and Ub are free pa- 
rameters in our model, and £,i,,o denotes the thermal energy density 
at the cluster centre. Rather than applying a scaling with the gas 
density as non-radiative MHD simulations by Dolag et al. (1999, 
2001) suggest, we chose the energy density of the thermal gas. This 
quantity is well behaved in the centres of clusters where current 
cosmological radiative simulations, that do not include radio-mode 
feedback from AGN, have an over-cooling problem which results 
in an overproduction of the amount of stars, enhanced central gas 
densities, too small central temperatures, and too strong central en- 
tropy plateaus compared to X-ray observations. Theoretically, the 
growth of magnetic field strength is determined through turbulent 
dynamo processes that will saturate on a level which is determined 
by the strength of the magnetic back-reaction (e.g., Subramanian 
2003; Schekochihin & Cowley 2006) and is typically a fraction of 
the turbulent energy density that itself should be related to the ther- 
mal energy density, thus motivating our model theoretically. 



3 RESULTS 

3.1 Radio synchrotron emission 

3.1.1 Simulated synclirotron scaling relations 

In order to determine the cluster scaling relations for non-thermal 
observables, we integrate the total surface brightness (composed 
of primary and secondary emission components) within the virial 
radius of each cluster. In our radiative simulations, we cut the 
region with r < 0.025 i^vir around the brightest central point- 
source that is caused by over-cooling gas of the cD galaxy. Since 
the modelled non-thermal emission processes reflect active non- 
equilibrium structure formation processes, we expect a large scat- 
ter in these scaling relations. Ideally, we would like to have a large 
sample of independent clusters to obtain reliable measurements of 
the scaling parameters. Thus, our limited sample will have larger 
uncertainties in the derived parameters. Figure I shows our simu- 
lated synchrotron scaling relations at y = 1.4 GHz using the total 
radio synchrotron luminosity within R^\y of all clusters. We note that 
the radio emission volume is significantly enlarged for our merg- 
ing clusters mostly due to the larger contribution of primary radio 
emission in the cluster outskirts (Paper II). To simplify comparison 
with observed giant radio halo samples, we additionally fit radio 
synchrotron scaling relations for our subsample of eight merging 
clusters. The fit parameters for our models with varying simulated 
physics and magnetic parameters can be found in Table 3. The fol- 
lowing conclusions can be drawn. 

Contributions of different emission components. (I) The 
secondary emission component is dominant for relaxing cool core 
(CC) clusters, and those that only experience a minor merger. The 
primary component exceeds the secondary one for major merging 
clusters by a factor of four as can be seen in our large post-merging 
cluster g72a with M =^ lO'^A"' Mg. (2) The secondary radio emis- 
sion is remarkably similar for our massive clusters while the scatter 
of the secondary emission increases notably for our small clusters 
with M < 2 X 10'"* Mq. This is due to the property of the hier- 
archical scenario of cluster formation which implies that virtually 
every large cluster is formed through a series of mergers of smaller 
progenitors. Each of these merging events triggered violent shock 
waves that accelerated CR protons through diffusive shock accel- 
eration. Over its cosmic history, these CRs accumulated within the 




51, B„ = lOuC aj = 0.5 

52, B„= lOnG, aj = 0.5 
S2, B„= 3nG, aj = 0.5 
S2, B„= lOuC aj = 0.25 

observed giant radio halos 
observed radio mini-halos 




in"'' ^ 

: + total synch., merger 

^ total synch., non-merger 

10"^^ - O primary synch. 

: O secondary synch. 



X 10-11 ^ 




10" 



M,„[;r' Ms] 



Figure 1. Cluster scaling relations for the radio synchrotron luminosities 
at y = 1 .4 GHz. The top panel shows the scatter of the individual clusters 
for our non-radiative (model SI) and radiative simulations (model S2). The 
middle panel shows the dependence of the scaling relations on the uncer- 
tainty in the magnetic field and simulated physics. The bottom panel shows 
the contribution of the individual emission components (primary, secondary 
radio synchrotron emission) to the total radio luminosities in our model S2 
while assuming a central magnetic field strength of So = 10 and an 
energy density scaling of ag = 0.5. 



cluster volume due to their cooling time being longer than the Hub- 
ble time (Volk et al. 1996; Berezinsky et al. 1997). The secondary 
radio emission probes the CR proton pressure which traces the time 
integrated non-equilibrium activities of clusters and is only modu- 
lated by the recent dynamical activities (see also Paper I, for aver- 
age values of the relative CR energy in different dynamical cluster 
environments). In our less massive clusters, the larger scatter of the 
secondary emission level is due to the larger variation of merging 
histories of these clusters and their weaker gravitational potential. 
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Table 3. Cluster scaling RELAnoNS for non-thermal observables*''. 

radio synchrotron: IC (£ic > 10 keV): y-rays {E^ > 100 MeV): 

aU clusters merging clusters 



model<2) 


g(3) 


"b 


v,0 


j8v 


v,0 




AS) 

-^c,o 


Ac 


/■(6) 

y,o 




SlP) 


10 


0.5 


0.78 ± 0.05 


1.45 ±0.07 


0.80 ±0.12 


1.50 ±0.17 


2.95 ± 0.26 


1.52 ±0.10 


7.85 ± 0.58 


1.52 ±0.09 




10 


0.5 


0.76 ±0.12 


1.32 ±0.18 


1.27 ± 0.26 


1.50 ±0.23 


1.66 ±0.22 


1.34 ±0.16 


5.46 ± 0.74 


1.32 ±0.16 


S2P) 


10 


0.25 


1.36 ±0.19 


1.30 ±0.16 


2.59 ± 0.41 


1.51 ±0.18 










SlP) 


3 


0.5 


0.14 ±0.01 


1.40 ±0.09 


0.26 ± 0.03 


1.46 ±0.20 










S2P) 


3 


0.5 


0.14 ±0.03 


1.22 ±0.23 


0.14 ±0.07 


1.46 ±0.30 


2.27 ± 0.29 


1.33 ±0.15 


5.65 ± 0.76 


1.32 ±0.16 


S2(2) 


3 


0.25 


0.31 ± 0.05 


1.27 ±0.18 


0.59 ±0.12 


1.49 ± 0.21 










S3P) 


10 


0.5 


0.89 ±0.16 


1.07 ± 0.21 


1.80 ±0.37 


1.31 ± 0.23 


2.24 ± 0.23 


1.24 ±0.12 


8.66 ± 1.12 


1.17 ±0.15 



Notes: 

(1) The cluster scaling relations for non-thermal observables are defined by A = AqM^j, where M15 = My„l{W^^'MQlh) and the respective non-thermal 
luminosity is obtained by integrating over the virial region of the cluster within if 200 and applying a central cut around the brightest central point-source for 
radii r < 0.025 Ryir- 

(2) The definition for our different models can be found in Table 2. 

(3) The definition for the parametrisation of the magnetic energy density is given by sb = eb,o (£th/£th,o)^"* according to (1) and Bq is given in units of |aG. 

(4) The normahsation of the radio synchrotron scaling relations is given in units of 10^^ erg s"' Hz"' hjo. 

(5) The normalisation of the IC scaling relations (Eic > 10 keV) is given in units of 10*' y s"' /170. 

(6) The normalisation of the y-ray scaling relations (Ey > 100 MeV) is given in units of 10*' y s"' /170. 



This leads to a larger modulation of the CR pressure and reflects 
inore sensitively the current merging activity of the cluster than it is 
the case in large systems. (3) In contrast to the secondary emission, 
the pressure of primary CR electrons sensitively resembles the cur- 
rent dynamical, non-equilibrium activity of forming structure and 
results in an inhomogeneous and aspherical spatial distribution with 
respect to collapsed objects. This leads to a large cluster-to-cluster 
variation of the primary radio emission. 

Normalisation: (1) The normalisation of the non-thermal 
scaling relations depends only weakly on whether radiative or non- 
radiative gas physics is simulated provided we consider in both 
cases only CRs from structure formation shocks. As discussed in 
Paper II, this is mainly due to self-regulated effects of the CR pres- 
sure. The CR cooling timescales due to Coulomb and hadronic in- 
teractions of CRs, Tpp/coui <^ igas. adjust to dilferent density levels in 
our simulations with radiative or non-radiative gas physics. Given 
a similar CR injection, this implies a higher CR number density for 
a smaller gas density hqr cc n^^.^^. The secondary CR emissivities 
(synchrotron, IC, or pion decay) scale as y'sec <^ ncR^gas <^ const and 
remain almost invariant with respect to different gas densities. (2) 
In contrast, the normalisation sensitively depends on the assump- 
tions and parametrisation of the magnetic field. This clearly shows 
the need to understand observationally how the properties of large 
scale cluster magnetic fields vary with cluster mass and dynamical 
state. 

Slope: (1) The slope of the radio synchrotron scaling relations 
for our merging cluster sample is largely independent of the sim- 
ulated physics or the parameters of our magnetic field if we only 
consider CRs from structure formation shocks (models S 1 and S2). 
The scaling relation is close to Lnt (details can be found 

in Tab. 3). The slope decreases to yS,, = 1.3 if we additionally ac- 
count for CRs from SNe feedback within galaxies. (2) If we con- 
sider all radio emitting clusters, i.e. we also account for radio-mini 
halos, the slope flattens in our radiative simulations by Afiy ^ 0.2. 
As a caveat, our scaling relations assume the same parametrisa- 
tion of the magnetic field for all clusters. If the central magnetic 
field scales with the cluster mass, the slopes will be accordingly 
steeper. Additionally, this self-similarity could be broken in the ra- 



dio synchrotron scaling relations, once inagnetic field are dynami- 
cally simulated and respond to the dynamical state of a cluster. 

Scatter: In our non-radiative simulations, the scatter in the 
radio synchrotron scaling relations is much smaller than in our ra- 
diative ones. There are no CC clusters in our non-radiative sim- 
ulations by definition. If a merger takes place, there are stronger 
shock waves in our radiative simulations due to the slightly cooler 
temperatures that imply smaller sound velocities and larger Mach 
numbers. Thus, the difi'erence between relaxed and merging cluster 
is more pronounced in our radiative simulations. 

3.1.2 Comparison to observations 

The observed sample of giant cluster radio halos (Cassano et al. 
2006) and that for cluster radio mini-halos (Gitti et al. 2004) is 
compared to our simulated scaling relations. 

Radio luminosity: Generally, our simulated giant cluster ra- 
dio halos show the same level of radio synchrotron emission as 
observed ones given a model of the magnetic field that is supported 
by Faraday rotation observations (Carilli & Taylor 2002; Widrow 
2002; Govoni & Feretti 2004, and references therein). 

Cluster magnetic fields: The radio synchrotron emissivity 
scales as 

i i-r c „(av+l)/2,,-ary 

7v'^£CRe«s V , (/) 

where ecrc and Eg denote the energy densities of CR electrons and 
magnetic fields, respectively, and the synchrotron spectral index 
ttv = aiiij/2 = (ffie - l)/2 is related to the spectral index of the in- 
jected electron population ttinj as well as to that of the cooled elec- 
tron population a^. Typical synchrotron spectral indices of clus- 
ter halos and relics span a range of = 1 . . . 1.3. This implies a 
similar contribution to the radio luminosity-mass scaling relation 
of clusters from the energy density of CR electrons and that of 
magnetic fields. Our radio synchrotron scaling relations assume the 
same physical model for the magnetic field irrespective of cluster 
mass and dynamical state. Conversely, we can interpret our simu- 
lated synchrotron scaling relations as tracks in the radio luminosity- 
cluster mass plane which are labelled with a set of parameters of 
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our magnetic model such as central magnetic field and magnetic 
decline. Radio mini-halos tend to have a higher radio luminosity 
compared to the giant radio halos. This hints towards a larger cen- 
tral magnetic field of the order of 10 \iG in these relaxed cool core 
clusters compared to the apparently preferred weaker central field 
strength of the order of 3 |J.G in merging systems. Interestingly, 
this characteristics of cluster magnetic fields is also consistent with 
Faraday rotation measurements (Vogt & EnBlin 2005, and refer- 
ences therein). Radio mini halos have been rarely observed in re- 
laxed cool core clusters. This might be partly due to the strong radio 
emitting AGN at the centres of cool core clusters which implies a 
large dynamic flux range to the underlying diffuse radio mini-halo 
and makes them very challenging to observe. 

Correlation between radio halos and mergers: How do our 
simulations support the observed radio halo-merger correlation? (1) 
The radio emission from primary, shock-accelerated electrons can 
boost the total radio emission of a major merging clusters by a fac- 
tor of four (cf. our massive post-merging cluster g72a). This factor 
sensitively depends on the mass ratio, geometry, and the advanced 
state of the merger. (2) In cool core as well as in merging clusters, 
the central radio emission is dominated by synchrotron emission 
from hadronically generated electrons. In Paper II, we show that 
the emission size of the simulated giant radio halos is increased 
due to the complex network of virializing shock waves in the clus- 
ter periphery that are able to efficiently accelerate CR electrons and 
amplify the magnetic fields due to strong shear motions. (3) The 
relative CR pressure is modulated by current merger activity of a 
cluster. While this positive modulation is small in massive systems, 
it can be substantial for less massive systems as can be inferred 
from Fig. I. This is due to larger variation of merging histories and 
the smaller gravitational potential in small clusters that causes the 
radio emission to respond more sensitively to merging activity. 

Observed scatter: The merger causes clusters that are hosting 
a giant radio halo to depart from hydrostatic equilibrium and leads 
to a complicated morphology that in general is not spherical. As a 
result, the masses in merging clusters can be either over- or under- 
estimated, depending on the amount of turbulent pressure support, 
the presence of shocks, and the amount of substructure which tends 
to flatten the average density profile (Evrard et al. 1996; Roettiger 
et al. 1996; Schindler 2002). The sample of observed giant radio 
halos of Cassano et al. (2006) scatters by two orders of magnitude 
in synchrotron luminosity while the virial masses of the hosting 
clusters only spans a factor of three. Barring observational uncer- 
tainties, the large range of d3mamical states and merger geometries 
among clusters as well as the variation of the magnetic proper- 
ties such as central field strength and magnetic decline furthermore 
contribute to the scatter in the scaling relations. The small sample 
size in combination with the mentioned uncertainties make it im- 
possible to determine a reliable observational synchrotron scaling 
relation for radio halos. Thus, the simulated scaling relations can 
only be compared to the total luminosity of the observed clusters. 
Studies of radio synchrotron emission from clusters should be com- 
plemented by studies of the pixel-to-pixel correlation of the syn- 
chrotron and X-ray surface brightness (Govoni et al. 2001; Pfrom- 
mer et al. 2007b). 

3.2 Inverse Compton and pion decay induced 7-ray emission 

In contrast to the observed diffuse radio synchrotron emission from 
clusters, y-rays from clusters have not been detected yet (Reimer 
et al. 2003). In principle, inverse Compton and pion decay induced 
y-ray emission are the cleanest way of probing structure formation 



shock waves and the accelerated CR electron and proton popula- 
tions since these non-thermal emission processes are not weighted 
with the magnetic energy density as it is the case for synchrotron 
emission. 

3.2.1 Inverse Compton and y-ray cluster scaling relations 

We determine the cluster scaling relations for the non-thermal y-ray 
luminosities as before in Sect. 3.1.1. Figure 2 shows our simulated 
IC and pion decay scaling relations. The fit parameters for our mod- 
els with varying simulated physics and magnetic parameters can be 
found in Table 3. 

Contributions of different emission components: In the IC 
scaling relations (Eic > 10 keV), we see a similar picture as we 
found for the radio synchrotron scaling relations, albeit somewhat 
amplified since the weighting with the magnetic energy density is 
negligible at these energy bands. The secondary emission com- 
ponent is dominant for relaxing CC clusters, and those, that only 
experience a minor merger. The primary component exceeds the 
secondary one for major merging clusters. In the y-ray scaling re- 
lations (Ey > 100 MeV), the pion decay component is always 
dominant over the primary and secondary IC emission components. 
This finding does only weakly depend on the assumed spectral in- 
dex for the CR proton distribution function since the energy band 
Ey > 100 MeV is dominated by the peak of the pion bump that is 
produced by GeV-protons (Pfrommer & EnBlin 2004a). 

Normalisation: (1) The normalisation of the non-thermal 
scaling relations depends only weakly on whether radiative or non- 
radiative gas physics is simulated provided we consider in both 
cases only CRs from structure formation shocks. As previously dis- 
cussed (cf. Sect. 3.1.1), this is mainly due to self-regulated effects 
of the CR pressure due to CR cooling mechanisms. (2) If we ad- 
ditionally account for CRs from SNe feedback within galaxies, the 
normalisation increases due to the second source of CR injection. 
This increase is higher for our y-ray scaling relations which are 
completely dominated by the pion decay emission component. To 
which extend CRs are able to diffuse out of the cold ISM and enrich 
the ICM needs to be studied separately. 

Slope: The slope of the non-thermal IC/y-ray scaling rela- 
tions depends weakly on the simulated physics and is almost in- 
dependent of the parameters of our magnetic field. For all three 
non-thermal emission mechanisms (synchrotron, IC, pion decay 
induced y-ray emission), very similar slopes are found. This is a 
non-trivial finding, since the relative contribution of the various 
emission components differs for the different energy bands consid- 
ered in this paper. Our set of non-radiative simulations (SI) yields a 
slope of j8ic,y - 1.5. This reduced in our radiative simulations (S2) 
to jSicy - 1.33 and furthermore decreased when considering CRs 
from SNe feedback (S3) to /3ic.y ^ 1.2. 

Scatter: In our non-radiative simulations, the scatter in the y- 
ray scaling relations is somewhat smaller than in our radiative ones 
while it is similar in the IC scaling relations. There are two rea- 
sons for this. (1) In our non-radiative simulations, there are no CC 
clusters by definition. In merging clusters, there are stronger shock 
waves in our radiative simulations due to the slightly cooler tem- 
peratures that imply smaller soimd velocities and larger Mach nimi- 
bers. This leads to more effective diflFusive shock acceleration and 
an enhanced level of non-thermal emission. (2) The primary emis- 
sion component has its largest impact for the IC hard X-ray emis- 
sion (compared to the y-ray emission). This component is largely 
responsible for the large scatter since it traces the current dynami- 
cal, non-equilibrium activity of the cluster. Looking at the individ- 
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IC emission (Ey > 10 keV): j-ray emission (Ey > 100 MeV): 




Figure 2. Cluster scaling relations for non-thermal observables. Shown are the relations for the inverse Compton (IC) hard X-ray luminosities for Ey > 10 keV 
(left panel) and the relations for the y-ray luminosities for Ey > 100 MeV (right panel). The top panels show the scatter of the individual clusters for our 
non-radiative (model SI) and radiative simulations (model S2). The middle panels show the dependence of the scaling relations on the uncertainty in the 
magnetic field and simulated physics. The bottom panels show the contribution of the individual emission components (primary, secondary, pion decay y-rays) 
to the total cluster luminosities in our model S2 while assuming a central magnetic field strength of Bo = 10 M^G and an energy density scaHng of ag = 0.5. 



ual non-thermal luminosities of our clusters (top panels in Fig. 2), 
one can notice a large scatter. In particular for the y-ray emission, 
this scatter increases for less massive clusters in our radiative mod- 
els and can boost the y-ray luminosity up to a factor of four. Due to 
the small sample size of our simulated high-resolution clusters, we 
are unable to statistically quantify this effect reliably. 



3.2.2 Luminosity and flux functions 

We combine our derived cluster scaling relations for non-thermal 
observables with the complete sample of the X-ray brightest galaxy 
clusters (HIFLUGCS, Reiprich & Bohringer (2002)) to predict 
fluxes and luminosities of each of the clusters. For the luminos- 
ity distance and the cluster masses, we assumed a ACDM cos- 
mology with a currently favoured Hubble constant, hjQ, where 
Hq = lOhjo km s"' Mpc"'. Dependent on the simulated physics, 
we can thus derive flux and luminosity functions for the y-ray emis- 
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Figure 3. We use the complete sample of the X-ray brightest clusters (HIFLUGCS, Reiprich & Bohringer (2002)) to predict flux and luminosity functions of 
the y-ray emission for Ey > 100 MeV. The definition for our different models can be found in Table 2. The top panels shows the differential flux/luminosity 
functions while the bottom panels show the respective cumulative functions. Assuming a GLAST sensitivity after two years of 2 X 10~'7 cm~^ s"' , we predict 
the detection of seven to eleven clusters named in the top left panel, depending on the adopted model. 



sion (Fig. 3) and for the hard X-ray IC emission (Fig. 4). Tables Al 
and A2 in the Appendix show the 15 brightest as well as the 15 
most luminous clusters of our homogeneous flux limited sample. 

7-ray emission: Assuming a GLAST sensitivity after two 
years of 2 x lO^'y cm"^ s"', we predict the detection of seven to 
eleven clusters named in the top left panel of Fig. 3, depending on 
the adopted model. The brightest y-ray clusters are Ophiuchus, For- 
nax, Coma, A3627, Perseus, and Centaurus (A3526), independent 
of the simulated physics. Among these, only Ophiuchus is among 
the ten most massive and thus most luminous clusters of the HI- 
FLUGCS sample. This statement can be modified once we consider 
scatter in y-ray luminosity due the the varying dynamical states of 
these clusters which might modify the rank ordering of the individ- 
ual systems. 

The distribution of the number of clusters with a given y-ray 
flux Ty is flat in the variable logTy down to '7^ 10^'y cm"^ s"' 
where the true number of clusters is suddenly increasing towards 
fainter fluxes. The large intrinsic scatter around the scaling rela- 
tion, especially at low y-ray luminosities, is expected to increase 
the number of cluster detections for GLAST by scattering clusters 
above the survey flux limit for our case of a decreasing differential 
distribution with increasing flux.' For comparison, we show how 



This effect is also known as Eddington bias (Eddington 1913). 



a flat differential distribution with a maximum flux !Fmax translates 
into a cumulative one (dotted line): 

N(>r) = No f 0(log T^,, -x)dx = No log ^ , (3) 

Jiogr 

where 9(x) denotes the Heaviside function. Any cumulative flux 
function steeper than the dotted line benefits from the scatter 
around the scaling relation. This has to be taken into account when 
deriving the observed y-ray luminosity function. However, due 
to the limited statistics in our simulated sample, further work is 
needed to quantify this expected scatter. 

The luminosity function shows an exponential cutoff at high 
y-ray luminosities that is inherited from the Press-Schechter mass 
function. The uncertainty at the high-mass end of our scaling rela- 
tions of our different CR models translates into a similar uncertainty 
of the exponential cut-off of the cumulative luminosity function. 
The decrease of the luminosity function at smaller luminosities is 
due to the incompleteness of the X-ray flux limited cluster sample. 

Inverse Compton emission: We predict Ophiuchus to be the 
brightest hard X-ray emitting cluster with a photon number flux of 
(1.7 . . .3.4)x 10^'y cm"^ s"' for energies > 10 keV. Despite the 
fact that the derived slope jSic of the IC scaling relation of model S3 
is different compared to the y-ray case, all our brightest y-ray clus- 
ters remain the brightest hard X-ray emitting clusters. The slopes 
in our other models are identical between the y-ray and IC case 
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IC hard X-ray flux function (Ey > 10 keV): 



IC luminosity function (Ey > 10 keV): 




'Fy [y cm ^ s 



Figure 4. We use the complete sample of the X-ray brightest clusters (HIFLUGCS, Reiprich & Bohringer (2002)) to predict flux and luminosity functions 
of the hard X-ray IC emission for Ey > 10 keV. The definition for our different models can be found in Table 2. The top panels shows the differential 
flux/luminosity functions while the bottom panels show the respective cumulative functions. 



which leads to identical rank ordering of the IC emitting clusters. 
Our findings with respect to the selection bias of the number of de- 
tected clusters due to the scattering in the scaling relations applies 
also in this case. 



3.2.3 Comparison to obserx'ations and previous work 

7-ray emission: Reimer et al. (2003) derived the EGRET upper 
limits on the high-energy y-ray emission of galaxy clusters us- 
ing nine years of successive observations. Stacking a sample of 
58 clusters and carefully accounting for the diffuse y-ray back- 
ground yielded an upper 2cr limit for the average cluster of 6 x 
lO^'y cm"- s"' for Ey > 100 MeV. The limits on individual clus- 
ters that this work predicts to be the brightest y-ray emitters range 
within (3 . . . 5)x lO^^y cm"^ s"' . Our predicted fluxes are consistent 
with these upper limits, providing an important consistency check 
of our models. 

Inverse Compton emission: There seems to be growing evidence 
for an excess of hard X-ray emission compared to the expected 
thermal bremsstrahlung in a number of clusters that is based on 
observations with instruments on board five different X-ray satel- 
lites. Prominent examples include the Coma cluster (Rephaeli et al. 
1999; Fusco-Femiano et al. 1999; Rephaeli & Gruber 2002; Fusco- 



Femiano et al. 2004, 2007b;^ Eckert et al. 2007; using the Rossi X- 
ray Timing Explorer (RXTE), BeppoSAX, and INTEGRAL) and the 
Perseus cluster (Sanders et al. 2005; Molendi 2007, using Chan- 
dra and XMM-Newton). Using our simulations, we can test the 
currently favoured hypothesis that this emission is due to inverse 
Compton radiation by CR electrons. Fusco-Femiano et al. (1999) 
claimed an excess flux of 2 x 10"" erg cm"' s"' between Ei = 
20 keV and Ei = 80 keV. For the Coma cluster, our models predict 
a inverse Compton number flux of ( 1 .3 ... 2.3) x 10"''7 cm"- s"' for 
energies above iiic,o =10 keV. To relate the number flux to an en- 
ergy flux, we assume a photon index of o-y = 1.15 and a scaling 
of Tic = To (£/£ic,o)"''- Using the notation for energy and number 
fluxes described in Paper II, we can calculate the energy flux in the 
observational hard X-ray band. 



E2 

ElC,Q 



(4) 



4 X 10 erg cm ^ s '. 



This is a factor of 50 below the claimed detection of hard X-ray 



^ The results of these papers have been challenged by an analysis that takes 
into account all systematic uncertainties in the critical pai'ameters including 
the choice of a source-free background field and the modelling of the ther- 
mal model for the ICM (Rossetti & Molendi 2004, 2007). 
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emission. We will discuss the implications of this discrepancy be- 
low. 

The same argument applies to the hard X-ray emission in 
the Perseus cluster where Sanders et al. (2005) find a flux of 
6.3 X 10 " erg cm"^ s ' between 2 and 10 keV. Assuming a pho- 
ton index of Uy - 1.15, this flux exceeds our IC prediction of 
5x 10"'^ erg cm~^ s"' for the same energy range by over two orders 
of magnitudes. We note that in the particular case of Perseus, the 
main cluster temperature of kT^ = 7 keV (Churazov et al. 2003) is 
very close to the energy limit of Chandra's imaging spectrometer, 
leaving a small lever arm for the detection of power-law compo- 
nent on top of the expected thermal bremsstrahlung components. 
Assuming that the hard X-ray emission is due to IC emission of 
CR electrons, we expect the non-thermal emission to be physically 
and spatially unrelated to the thermal emission components. The 
morphological similarity of the high energy/temperature emission 
maps of Sanders & Fabian (2007) show a clear spatial (or angular) 
anti-correlation between the hottest thermal components (4, 8 keV) 
and the power-law component which questions the IC interpreta- 
tion of the data. A thorough covariance analysis of the seven differ- 
ent emission components in the model of Sanders & Fabian (2007) 
would be needed in order to exclude the possibility of component 
coupling in their spectral deconvolution procedure that mimics a 
non-present power-law component. The IC interpretation is also 
challenged on theoretical grounds since it requires the energy den- 
sity of CR electrons to be in equipartition with the thermal plasma, 
leaving no room for rclativistic protons that have a much longer 
lifetime compared to electrons. 

Previous work: Most of the previous work that calculated the y- 

ray emission from individual clusters made very simplifying as- 
sumptions about the amount and spatial distribution of CRs within 
galaxy clusters (for a comprehensive review, cf. Blasi et al. 2007). 
Based on simplified analytical arguments such as spherical geome- 
try, virial equilibrium, and CRs that are diffusing from a source in 
the cluster centre, Colafrancesco & Blasi (1998) derive a scaling 
relation of the hadronically induced y-ray luminosity with cluster 
mass Jloc M^l^ that is much shallower than our relations in Table 3. 
The difference can be easily explained by our more realistic simula- 
tions that self-consistently follow the relevant CR physics leading 
to an inhomogencous distribution of rclativistic protons, include 
hydrodynamical non-equilibrium effects and arbitrary cluster ge- 
ometries, and account for reaUstic cosmological merger histories. 

Modelling the non-thermal emission from clusters by numer- 
ically modelling discretised CR energy spectra on top of Eulerian 
grid-based cosmological simulations, Miniati et al. (2001a,b) de- 
rive various scaling relations of non-thermal cluster emission rang- 
ing from radio synchrotron, IC soft and hard X-rays, to y-rays 
which are in part considerably steeper than our relations in Ta- 
ble 3. In contrast to our approach, these models neglected the hy- 
drodynamic pressure of the CR component, were quite limited in 
their adaptive resolution capability, and they neglected dissipative 
gas physics including radiative cooling, star formation, and super- 
nova feedback. The cluster sample was comprised of small systems 
with average core temperatures of 0.3 keV < kT < 3 keV and non- 
thermal luminosities have been computed within a fixed radius that 
various between 1.5i?vii and 4i?vir for the smallest groups where 
/f2oo - 300 kpc. The discrepancy of the non-thermal scaling re- 
lations can be understood by two main effects that lead to an overes- 
timation of the CR pressure inside the clusters simulated by Miniati 
et al. (2001a) and thus overproduced the resulting non-thermal 
emission particularly in larger systems: (1) Miniati et al. (2000) 
identified shocks with Mach numbers in the range 4 < Al < 5 



as the most important in thermalizing the plasma. In contrast, Ryu 
et al. (2003) and Pfrommer et al. (2006) found that the Mach num- 
ber distribution peaks in the range 1 < < 3. Since diffusive 
shock acceleration of CRs depends sensitively on the Mach num- 
ber, this implies a more efficient CR injection in the simulations 
by Miniati et al. (2001a). (2) The grid-based cosmological simu- 
lations have been performed in a cosmological box of side-length 
50/!"' Mpc with a spatial resolution of 200/;"' kpc, assuming an 
Einstein-de Sitter cosmological model (Miniati et al. 2001a). The 
lack of resolution in the observationally accessible, dense central 
regions of clusters in the grid-based approach underestimates CR 
cooling processes such as Coulomb and hadronic losses. Secondly, 
these simulations are unable to resolve the adiabatic compression 
of a composite of CRs and thermal gas, an effect that disfavours the 
CR pressure relative to the thermal pressure. 



3.3 Minimum 7-ray flux 

For clusters that host giant radio halos with an observed luminos- 
ity Ly, we are able to derive a minimum 7-ray flux in the hadronic 
model. The non-detection of y-ray flux below this flux level lim- 
its the contribution of secondary radio emission to the giant radio 
halo independent of the spatial distribution of CRs and thermal gas. 
The idea is based on the fact that the radio luminosity of an equi- 
librium distribution of CR electrons, where injection and cooling 
is balanced, becomes independent of the magnetic field in the syn- 
chrotron dominated emission regime for sb » Cph (cf. Fig. 3 in 
Paper II), 



\(av-l)/2 



Ly = 



Ay J dVCpHN, for Sb » Sph and Oy ~ 1, 



dVCpMN, 



(5) 

(6) 



where Ay and Ay are constants of the hadronic interaction physics 
and given in the Appendix of Paper II, the volume integral ex- 
tends over the entire cluster, Cp oc m^r is the normaUsation of 
the CR momentum distribution and proportional to the CR num- 
ber density, is the number density of target nucleons for the 
hadronic interaction, gph = Ecmb + Cstars is the energy density of 
the cosmic microwave background (CMB) and the starlight pho- 
ton field, where the equivalent magnetic field strength of the en- 
ergy density of the CMB is given by Scmb = 3.24 ^i.G (1 + z)^, and 
Sbc - 31 (v/GHz) |i.G is a frequency dependent characteristic mag- 
netic field strength for synchrotron radiation. In this strong field 
limit, the volume integral of the synchrotron emission is equal to 
that of the y-ray emission resulting from pion-decay and can be 
eliminated yielding 



(7) 



and Z)iu,^ is the luminosity distance to the cluster. Smaller mag- 
netic fields would require a larger energy density of CR electrons 
in order to reproduce the observed synchrotron emission and thus 
enhance the simultaneously produced y-ray emission. For the sam- 
ple of known giant radio halos (Cassano et al. 2006), the Coma 
cluster is expected to have the largest y-ray flux since the combina- 
tion Lyobs/i'funi 1^^^' f'^™' times larger than that in other cluster 
that are hosting giant radio halos. The lowest possible hadronic y- 
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ray flux is realised for iiard CR spectral indices, ttp = 2, yielding 
^,min = 7.5 X 10""7 cm"^ s"' in Coma. 

It turns out, that this limit can be considerably tightened by 
requiring the average magnetic energy density to be locally less 
than the thermal energy density. For our Coma-like cluster g72a in 
our simulation, a central magnetic field strength of 10 (iG corre- 
sponds to a ratio of thermal-to-magnetic pressure of 20. Since the 
thermal pressure decreases by two orders of magnitude towards the 
virial radius, a constant magnetic energy density (as required by the 
synchrotron dominated emission regime) would exceed the thermal 
energy density by a factor of five. This requires knowledge of the 
spatial distribution of CRs, magnetic fields and thermal gas in our 
Coma-Uke cluster simulation, 

-T y,mm — . A r\2 ' -' 

where Ly,g72a is the central radio halo emission due to hadroni- 
cally produced CR electrons in our model S2 (CR acceleration 
at structure formation shocks while allowing for all CR loss pro- 
cesses). The predicted y-ray luminosity in this model amounts to 
-Cy,g72a - 7.3 X lO"" y s"' and is weakly dependent on the as- 
sumed CR spectral index of a = 2.3. The very conservative y- 
ray limit assumes a central magnetic field Bq = 10 H-G, ensures 
Pth > 2Pb everywhere within the virial region of the cluster and 
yields T'y,mm = 4x 10"'°y cm"^ s ' . For the same Bo and Pa, > 20^^ 
at the virial radius, we obtain If^min = 9 x lO'^^y cm~^ s"'. Adopt- 
ing an even lower central magnetic field fio - 3 ]xG as Faraday 
rotation studies of the Coma cluster indicate (Kim et al. 1990) 
and requiring Pa, > 20Pb at the virial radius, we obtain T^^mm = 
2x lO^'y cmr^ s"' = Tglast, 2yr, i.e. the GLAST all-sky survey will 
be able to scrutinise this scenario after two years. We would like to 
close this section by noting that our simulations predict a y-ray 
flux from Coma of Ty,mia = (4 ... 7) x 10"'y cm~^ s"'. This in turn 
would imply a central magnetic field Bq ^ 3 \iG with a constant 
average ratio of thermal-to-magnetic pressure of 200, comparing to 
the observed synchrotron flux. 



4 DISCUSSION AND CONCLUSIONS 

We performed high-resolution simulations of a sample of 14 galaxy 
clusters that span a mass range of almost two orders of magnitude 
and follow self-consistent CR physics on top of the dissipative gas 
physics including radiative cooling, star formation, and supernova 
feedback. The modelled CR physics in our simulations and our on- 
the-fly identification scheme of the strength of structure formation 
shock waves allows us to reliably compute the relativistic elec- 
tron populations at high energies. We consider relativistic electrons 
that are accelerated at cosmological structure formation shocks (so- 
called primary electrons) and those that are produced in hadronic 
interactions of cosmic rays with ambient gas protons (hence the 
name secondary or hadronic electrons). 



4.1 Non-thermal scaling relations 

In this paper, we concentrate on three observationally motivated 
wave-bands. (1) Radio synchrotron emission at 1.4 GHz, (2) non- 
thermal hard X-ray emission at energies Ey > 10 keV, and (3) y-ray 
emission at energies Ey > 100 MeV. We study the contribution of 
the different emission components to the total cluster luminosity in 
each of these bands, derive cluster scaling relations, and study their 



dependence on the simulated physics and adopted parametrisation 
of the magnetic field. Our main findings are as follows: 

(i) The secondary emission component (radio synchrotron and 
inverse Compton) is dominant for relaxing cool core clusters, and 
those that only experience a minor merger. The primary component 
can exceed the secondary one for major merging clusters by a factor 
of four. In the y-ray scaling relations (Ey > 100 MeV), the pion de- 
cay component is always dominant over the primary and secondary 
IC emission components. 

(ii) The normalisation of the non-thermal scaling relations de- 
pends only weakly on whether radiative or non-radiative gas 
physics is simulated provided we consider in both cases only 
CRs from structure formation shocks. This is mainly due to self- 
regulated effects of the CR pressure due to the density dependent 
CR cooling mechanisms. In contrast, the normalisation of the ra- 
dio synchrotron scaling relation sensitively depends on the assump- 
tions and parametrisation of the magnetic field. This clearly rein- 
forces the need to understand observationally how the properties 
of large scale cluster magnetic fields vary with cluster mass and 
dynamical state. 

(iii) The slope of the non-thermal scaling relations depends 
weakly on the simulated physics and is almost independent of the 
parameters of our magnetic field. For all three non-thermal emis- 
sion mechanisms (synchrotron, IC, pion decay induced y-ray emis- 
sion), very similar slopes are found. This is a non-trivial finding, 
since the relative contribution of the various emission components 
differs for the different energy bands considered in this paper. Our 
set of non-radiative simulations (SI) yields a slope of fSjcy - 1-5. 
This is reduced in our radiative simulations (S2) to /3ic,y - 1.33 
and furthermore decreased when considering CRs from SNe feed- 
back (S3) to jSic.y - 1.2. The slope of the synchrotron scaling re- 
lation steepens if we only consider merging galaxy clusters. As a 
caveat for our synchrotron scaling relations, we assume the same 
parametrisation of the magnetic field for all clusters. If the cen- 
tral magnetic field scales with the cluster mass, the slopes will be 
accordingly steeper Additionally, this self-similarity could be bro- 
ken in the radio synchrotron scaling relations, once magnetic field 
are dynamically simulated and respond to the dynamical state of a 
cluster This reinforces the need to understand the observed scal- 
ing properties of the magnetic field in clusters before we can draw 
strong conclusions about the theory underlying the cluster radio 
emission. 

(iv) In our non-radiative simulations, we observe large scatter 
in all non-thermal scaling relations. This is mostly driven by active 
merging systems that trigger violent shock waves and thus boost the 
primary emission signal. Our results hint at a larger contribution of 
the scatter towards less massive systems due to their smaller grav- 
itational potential which needs to be checked with a larger cluster 
sample size. The large scatter will have important implications for 
the number of detectable y-ray emitting clusters by GLAST. 

4.2 Radio synchrotron emission 

The unified model of radio halos and relics has been put forward 
in our companion paper (Paper II) and is based on studies of the 
morphology, profiles, and expected polarisation of our simulated 
diffuse cluster radio synchrotron emission. The derived radio lumi- 
nosities of the primary and secondary electron populations com- 
plement this picture. We are summarising the main findings of this 
work in the following. 

(i) Assuming magnetic field strengths provided by Faraday ro- 
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tation observations, we are able to successfully reproduce the ob- 
served radio synchrotron luminosities of giant radio halos as well 
as radio mini-halos in our simulations. 

(ii) Each of our radio halo scaling relations assumes one physi- 
cal model for the magnetic field that is described by a central field 
strength and magnetic decline. We assume it to be independent of 
the cluster mass and dynamical state. In this respect, the simulated 
scaling relations can be understood as contour lines in the radio 
luminosity-cluster mass plane which are labelled with a set of pa- 
rameters of our magnetic model. Radio mini-halos have a higher 
radio luminosity on average compared to that of giant radio halos. 
This points towards a larger central magnetic field of the order of 
10 \xG in these relaxed cool core clusters compared to the appar- 
ently preferred weaker central field strength of the order of 3 ^G in 
merging systems. This finding is consistent with Faraday rotation 
measurements and strongly hints at an amplifying mechanism of 
magnetic field strengths in relaxed clusters such as adiabatic com- 
pression of the fields during the formation of the cool core or AGN 
feedback amplified fields as argued in Enl51in & Vogt (2006). 

(iii) Observed giant radio halos are all associated with merging 
clusters. The merger causes these systems to depart from hydro- 
static equilibrium and leads to a complicated non-spherical mor- 
phology. The resulting X-ray mass estimates are subject to large 
uncertainties and might be partly responsible for the large scatter 
of observed giant radio halos of scatters by two orders of magni- 
tude in synchrotron luminosity while the virial masses of the host- 
ing clusters only spans a factor of three. Clustcr-to-cluster varia- 
tions of the geometry, mass ratio, the advanced state of the merger, 
and magnetic field strengths contribute furthermore to the scatter 
in the scaUng relations. The smaU sample size in combination with 
the mentioned uncertainties doom all attempts to determine a reli- 
able observational synchrotron scaling relation for radio halos. In 
contrast, studies of the pixel-to-pixel correlation of the synchrotron 
and X-ray surface brightness enable valuable insights that are not 
subject to the assimiption of spherical symmetry. 

4.3 Inverse Compton emission 

Given our reliable modelling of the synchrotron and IC emitting 
high-energy CR electron populations and our convenient parametri- 
sation of the magnetic field that is calibrated against Faraday ro- 
tation measurements, we can successfully reproduce the luminos- 
ity of observed giant radio halos. However, our predicted inverse 
Compton flux for the Coma and Perseus cluster falls short of the 
detected excess of hard X-ray emission compared to the expected 
thermal bremsstrahlung by at least a factor of 50. Lowering the 
magnetic field strength will not reconcile this discrepancy, since 
the IC emissivity of a steady state electron population is indepen- 
dent of the magnetic energy density in the low-field regime for 
B <K ficMB = 3.24 \iG. This finding can be rephrased as follows. 
Combining the observed diffuse radio synchrotron and IC emission 
allows to eliminate the ab initio unknown energy density of rela- 
tivistic electrons and to obtain an estimate for the magnetic field 
strength that typically reaches values of ~ 0.3 \iG (e.g., EnBlin & 
Biermann 1998; EnBlin et al. 1999). There are now three problems 
associated with these low field strengths that challenge the standard 
inverse Compton interpretation of the hard X-ray excess emission. 
(1) These field strengths are an order of magnitude smaller than 
those derived from Faraday rotation measurements which translates 
into two orders of magnitude in energy density. (2) The energy den- 
sity of CR electrons ScRe that is in turn needed to explain the radio 
halo emission would thus be two orders of magnitudes larger than 



what our model of the primary and secondary electron populations 
predict. (3) If we increased scRe by two orders of magnitude (due 
to a different injection mechanism such as AGN jets), the result- 
ing IC emission at Ey > 100 MeV would challenge upper limits 
on the y-ray emission imposed by EGRET (Reimer et al. 2003). 
Furthermore, the acceleration mechanism would have to single out 
fermions over protons in order not to violate the EGRET bounds 
by overproducing the simultaneously produced y-rays from pion- 
decay. 

There have been suggestions in the literature to circumvent the 
first problem: EnBlin (2004) suggests that one could in principle 
reconcile the observed discrepancy of magnetic field estimates, if 
there is a significant difference between volume and CRe weighted 
averages. This would require an inhomogeneous magnetic energy 
distribution, an inhomogeneous distribution of the CR electrons, 
and an anti-correlation between these two. These conditions could 
be produced by physical mechanisms which produce inhomoge- 
neous or intermittent magnetic fields and at the same time anti- 
correlate the CRe density with respect to the magnetic energy den- 
sity. As a very plausible mechanism, he suggests synchrotron cool- 
ing in inhomogeneous magnetic fields that provides naturally the 
required anti-correlation. The hadronic model in conjunction with 
peripheral shock acceleration that we studied in our series of sim- 
ulations would provide a CR electron injection rate which is not 
correlated with the magnetic field strength, as would be required 
by the above explanation of the discrepancy of magnetic field es- 
timates by the two methods would work. In contrast to this, in the 
re-acceleration model one would expect a strong positive correla- 
tion of CRe and magnetic field strength, since magnetic fields are 
essential for the CR electron acceleration. Petrosian (2001) allevi- 
ates the difficulties with the low magnetic field strengths in the IC 
model by taking into account effects of observational selection bias 
and evoking non-standard assumptions of a non-isotropic pitch an- 
gle distribution as well as spectral breaks in the energy distribution 
of the relativistic electrons. 

What are the model uncertainties of our simulations that might 
boost the energy density of CR electrons thus circumventing the 
second problem? (I) The scatter in the IC scaUng relations seems 
only to be able to account for another factor of two, albeit the small 
sample size of our simulated high-resolution clusters makes it im- 
possible to statistically quantify the scatter reliably. (2) Adopting 
central magnetic field strength B « Scmb will only increase the IC 
emissivity that is emitted by a steady state electron population by 
another factor of two compared to our low-field case of B ^ Scmb 
(cf. Fig. 3 in Paper II). (3) Are there any other sources that inject CR 
electrons homogeneously throughout the cluster volvune and resup- 
ply them on a time scale shorter than their radiative cooling time of 
T ^ 10** yrs? CR diffusion out of AGN and radio galaxies will not 
reproduce the required homogeneous distribution of CR electrons 
in order to explain radio halos. Secondly, diffusion will lead to a 
narrow, steep profile of the CR electron energy density with a max- 
imum radius of ^J{R^) = V6/cTcooi = 14 kpc, asstmring a large CR 
diffusivity of k = 10^' cm^/s and a combined IC/synchrotron cool- 
ing time of T = 10* yr that corresponds to IC emitting electrons at 
10 keV with a Lorentz factor of y - 3 x 10^ and assuming a mag- 
netic field of 8 |iG. This is an order of magnitude smaller than the 
emission radius of the Perseus radio-mini halo and is much smaller 
than the emission size of giant radio halos. The re-acceleration 
model might in principle have the correct properties to explain the 
spatial and spectral electron distribution. As laid out above, it faces 
however severe problems in reconciUng the observed discrepancy 
of magnetic field estimates from Faraday rotation measurements 
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on the one hand and combining synchrotron and IC measurements 
one the other hand. (4) Our model of the diffusive shock acceler- 
ation mechanism assumes a featureless power-law (for details, see 
Paper II). Future work will be dedicated on improving this model to 
incorporate more elaborate plasma physical models. (5) In the liter- 
ature, the excess of hard X-ray emission compared to the expected 
thermal bremsstrahlung in the Coma cluster is currently controver- 
sially discussed (Fusco-Femiano et al. 2004, 2007a,b; Rossetti & 
Molendi 2004, 2007). Observational efforts, such as the future hard 
X-ray missions NuSTAR and Simbol-X, have to be undertaken to 
unambiguously detect the spectral and spatial characteristics of the 
hard X-ray excess emission. 

If one relaxes the requirement of explaining the hard X-ray 
excess with the same population of electrons that is responsible 
for the radio synchrotron emission, there are two other models that 
try to explain the non-thermal excess by synchrotron radiation of 
ultra-relativistic (multi-TcV) electrons and positrons. These elec- 
trons are continuously injected throughout the entire intra-cluster 
medium either by interactions of hypothetical very high-energy 
y-rays with diffuse extragalactic radiation fields (Timokhin et al. 

2004) or by means of pair production processes of CMB photons 
in the Coulomb field of ultra high-energetic CR (UHECR) pro- 
tons that are accelerated at structure formation shocks (Inoue et al. 

2005) . While the generation of extremely high-energy photons re- 
mains the main challenging question for the first model, the energy 
requirement of a UHECR population of 10*' erg is rather extreme. 
Both models, however, are not able to reproduce the radio-halo 
emission that is detected in these clusters. 

Possibly the most elegant explanation for the hard X-ray 
emission is non-thermal bremsstrahlung of a supra-thermal elec- 
tron population that is energised by Coulomb collisions be- 
tween the quasi-thermal electrons and non-thermal protons (Wolfe 
& Melia 2007). Such an electron population displays a higher 
bremsstrahlung radiative efficiency than a pure power law popu- 
lation thus avoiding the overheating problem of the thermal plasma 
(Petrosian 2001). The non-thermal protons would simultaneously 
be responsible for the Coma cluster's diffuse radio halo emission 
(within the unified scheme put forward in Paper II). 

4.4 High-energy 7-ray emission 

Our predicted y-ray fluxes of nearby galaxy clusters are consistent 
with EGRET upper limits of these clusters (Reimer et al. 2003). 
The brightest y-ray clusters are typically a factor of five smaller 
than the derived upper limits, which provides an important consis- 
tency check of our models. We note that our simulations have not 
been tuned to match these upper limits, instead we modelled the CR 
physics to our best knowledge and calculated the y-ray luminosity 
of our simulated clusters. 

We predict the detection of the pion decay induced y-ray emis- 
sion of seven to eleven galaxy clusters by GLAST, depending on 
the adopted model. The expected brightest y-ray clusters are Ophi- 
uchus, Fornax, Coma, A3627, Perseus, and Centaurus (A3526), in- 
dependent of the simulated physics. Due to the increasing slope of 
the differential cluster flux number distribution towards smaller y- 
ray fluxes and the large scatter in the scaling relations (especially 
for less massive systems), we expect the detected nmnber of clus- 
ters to increase somewhat since clusters are scattered above the sur- 
vey flux limit. For clusters that host giant radio halos, we are able to 
derive a minimum y-ray flux in the hadronic model independent of 
the spatial distribution of CRs and thermal gas. The radio luminos- 
ity of an equilibrium distribution of CR electrons, where injection 



due to hadronic CR interactions and cooling is balanced, becomes 
independent of the magnetic field in the S3'nclirotron dominated 
emission regime. A smaller magnetic field would require a larger 
energy density of CR electrons to reproduce the observed syn- 
chrotron luminosity and thus increase the associated y-ray flux. In 
Coma, the absolute minimum flux of fy^moi = 7.5x 10""y cm~^ s'^ 
is well below the sensitivity of GLAST. Assuming magnetic field 
strengths as derived by Faraday rotation measurements, these limits 
can be considerably tightened to match the GLAST sensitivity after 
two years of all-sky survey, IFglast, 2yr = 2 x 10^'y cm"^ s"'. The 
detection of hadronically induced y-ray emission will enable us to 
determine the CR proton pressure in clusters and unambiguously 
decide upon the model of cluster radio halos. 
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APPENDIX A: PREDICTING THE BRIGHTEST AND 
MOST-LUMINOUS NEARBY 7-RAY CLUSTERS 

High-energy y-ray fluxes and luminosities a related by the simple 
conversion formula. 



= 8.4x10-^2% f-_i2_)f ^'"^ \".(A1) 
cm2 s \ 10« s-i hio I \ 100 Mpc h.-^ j 

This paper has been t5?peset from a TgX/ 13T^ file prepared by the 
author. 
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Table Al. The brightest IC/y-RAY clusters of the HIFLUGCS sample:*'' 



cluster name 


z 


w(2) 
'"200 


^ic (SI) 


IC emission, 
^ic (S2) 


T7 ^ 1 n 

iLy > l\J KcV. 

T^^, (S2, B3) 


y^ic ^^^> 


T-y (SI) 


y-ray emission 
Ty (S2) 


y iLy > lUU iVlCV. 

Ty (S2, B3) 


Ty (S3) 


Ophiuchus 


0.0280 


2.32 


3.43 


1.77 


2.40 


2.26 


9.11 


5.75 


5.95 


8.49 


Fornax 


0.0046 


0.10 


1.12 


1.02 


1.42 


1.82 


3.04 


3.55 


3.68 


8.38 


Coma 


0.0232 


1.38 


2.30 


1.30 


1.78 


1.76 


6.12 


4.28 


4.43 


6.82 


A3627 


0.0163 


0.66 


1.51 


0.98 


1.35 


1.43 


4.04 


3.27 


3.39 


5.84 


Perseus 


0.0183 


0.77 


1.52 


0.96 


1.32 


1.38 


4.08 


3.20 


3.31 


5.57 


A3526 


0.0103 


0.27 


0.98 


0.75 


1.04 


1.21 


2.65 


2.56 


2.65 


5.21 


A 1060 


0.0114 


0.30 


0.96 


0.71 


0.99 


1.13 


2.57 


2.43 


2.51 


4.86 


M49 


0.0044 


0.05 


0.37 


0.39 


0.55 


0.76 


1.01 


1.38 


1.43 


3.67 


AWM7 


0.0172 


0.43 


0.72 


0.50 


0.69 


0.77 


1.94 


1.70 


1.76 


3.22 


3C129 


0.0223 


0.66 


0.81 


0.53 


0.72 


0.77 


2.18 


1.76 


1.82 


3.13 


NGC4636 


0.0037 


0.03 


0.20 


0.24 


0.34 


0.50 


0.56 


0.87 


0.90 


2.52 


A1367 


0.0216 


0.41 


0.41 


0.29 


0.40 


0.45 


1.10 


0.98 


1.02 


1.88 


A0754 


0.0528 


1.87 


0.67 


0.36 


0.49 


0.47 


1.78 


1.18 


1.22 


1.79 


Triangulum 


0.0510 


1.54 


0.54 


0.30 


0.41 


0.40 


1.43 


0.98 


1.01 


1.53 


NGC5846 


0.0061 


0.04 


0.14 


0.15 


0.22 


0.31 


0.39 


0.55 


0.57 


1.50 



Notes: 

(1) IC and y-ray fluxes of clusters that are contained in the complete sample of the X-ray brightest clusters (HIFLUGCS, the Highest X-ray FLUx Galaxy 
Cluster Sample, Reiprich & Bohringer (2002)). We predict these fluxes using our cluster scaHng relations for non-thermal observables defined in Table 3. The 
definition for our different models can be found in Table 2. The clusters are ordered according to their decreasing y-ray flux in oiu" model S3. 

(2) Mass contained within R200 m units of lO'^/i^Q M©. 

(3) Predicted total (primary and secondary) IC flux in units of 10"^y cm~^ If not otherwise mentioned, we use our magnetic parametrisation of 
Bo = 10 nG and as = 0.5. Our model with B3 refers to a smaller central value for the magnetic field of Bq = 3 |iG and yields a higher IC luminosity. 

(4) Predicted total (primary IC, secondary IC, pion decay) y-ray flux in units of lO'^ycm"^ '^^^ otherwise mentioned, we use our magnetic 
parametrisation of Bq = 10 ^G and as = 0.5. Our model with B3 refers to a smaller central value for the magnetic field of Bo = 3 (jG which barely effects the 
y-ray flux due to the dominant contribution from pion decay emission. 

Table A2. The most luminous IC/y-RAY clusters of the HIFLUGCS sample:*'' 











IC emission. 


£^ > 10 keV: 






y-ray emission, 


£^ > 100 MeV: 




cluster name 


z 


200 


X<c' (SI) 


Xic' (S2) 


X<c' (S2, B3) 


XJc' (S3) 


(SI) 


X^,"' (S2) 


X^f' (S2, B3) 


4^' (S3) 


A2163 


0.2010 


3.71 


12.64 


6.00 


8.10 


7.28 


3.35 


1.93 


1.99 


2.65 


A3888 


0.1510 


2.55 


7.14 


3.63 


4.92 


4.58 


1.89 


1.17 


1.21 


1.71 


A1914 


0.1712 


2.43 


6.62 


3.39 


4.60 


4.31 


1.76 


1.10 


1.14 


1.61 


Ophiuchus 


0.0280 


2.32 


6.16 


3.18 


4.32 


4.06 


1.64 


1.03 


1.07 


1.52 


A3827 


0.0980 


1.96 


4.77 


2.54 


3.46 


3.31 


1.27 


0.83 


0.86 


1.25 


A0754 


0.0528 


1.87 


4.45 


2.39 


3.25 


3.12 


1.18 


0.78 


0.81 


1.19 


A1689 


0.1840 


1.76 


4.06 


2.21 


3.00 


2.90 


1.08 


0.72 


0.75 


1.11 


A3266 


0.0594 


1.70 


3.83 


2.10 


2.85 


2.77 


1.02 


0.69 


0.71 


1.06 


A2065 


0.0721 


1.67 


3.74 


2.05 


2.79 


2.71 


0.99 


0.67 


0.69 


1.04 


A2256 


0.0601 


1.56 


3.36 


1.87 


2.54 


2.49 


0.90 


0.61 


0.63 


0.96 


Triangulum 


0.0510 


1.54 


3.30 


1.84 


2.50 


2.45 


0.88 


0.60 


0.62 


0.94 


A2142 


0.0899 


1.50 


3.18 


1.78 


2.42 


2.38 


0.85 


0.58 


0.60 


0.92 


A0644 


0.0704 


1.42 


2.91 


1.64 


2.24 


2.21 


0.78 


0.54 


0.56 


0.86 


Coma 


0.0232 


1.38 


2.81 


1.59 


2.17 


2.15 


0.75 


0.52 


0.54 


0.83 


A2029 


0.0767 


1.34 


2.68 


1.53 


2.08 


2.07 


0.71 


0.50 


0.52 


0.80 



Notes: 

(1) IC and y-ray luminosities of clusters that are contained in the complete sample of the X-ray brightest clusters (HIFLUGCS, the Highest X-ray FLUx 
Galaxy Cluster Sample, Reiprich & Bohringer (2002)). We predict these luminosities using our cluster scaling relations defined in Table 3. The definition for 
our different models can be found in Table 2. The clusters are ordered according to their decreasing y-ray luminosities in our model S3. 

(2) Mass contained within R200 in units of 10"/i:^q Mq. 

(3) Predicted total (primary and secondary) IC luminosity in units of lO'^'y s~'/j7o. If not otherwise mentioned, we use our magnetic parametrisation of 
Bo = 10 \iG and = 0.5. Our model with B3 refers to a smaller central value for the magnetic field of Bo = 3 nG and yields a higher IC luminosity. 

(4) Predicted total (primary IC, secondary IC, pion decay) y-ray luminosity in units of 10'*''y s^'Avq. If not otherwise mentioned, we use our magnetic 
parametrisation of Bq = 10 ^G and ag = 0.5. Our model with B3 refers to a smaller central value for the magnetic field of Bo = 3 [oG which barely efli'ects the 
y-ray luminosity due to the dominant contribution from pion decay emission. 
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